Pair Density Waves in coupled doped two-leg Ladders 



Javier Almeida, 1 Guillaume Roux, 2 and Didier Poilblanc 1 

1 Laboratoire de Physique Theorique, CNRS and Universite de Toulouse, UPS, (IRSAMC), F-31062 Toulouse, France 
2 Univ Paris-Sud, Laboratoire de Physique Theorique et Modeles Statistiques, 
UMR8626, Orsay F-91405, France; CNRS, Orsay, F-91405, France. 
(Dated: July 16, 2010) 

Motivated by Resonant X-ray scattering experiments in cuprate ladder materials showing charge order mod- 
ulation of period A = 3 and 5 at specific hole densities, we investigate models involving the electronic t-J 
ladders and bosonic chains coupled via screened Coulomb repulsion. Extensive density matrix renormalization 
group calculations applied to the ladders/chains supplemented by a self-consistent mean-field treatment of the 
inter-ladder/chain coupling provide quantitative estimates of the charge order for A = 3, 4 and 5. As previ- 
ously proposed, such patterns correspond to the emergence of pair density waves which stem from the strong 
electronic correlations. We comment on the existence of a A = 4 modulation not seen so far in experiment. 

PACS numbers: 71.45.Lr, 75.40.Mg, 74.72.Gh 



It is fascinating that electron correlations alone could lead 
to unconventional pairing, eventually leading to supercon- 
ducting or other exotic phases such as pair crystals, e.g. in 
the two-dimensional (2D) high critical-temperature cuprate 
superconductors. The quasi one-dimensional (ID) two-leg 
ladder material Sri^-rCa^O^O^ (SCCO) is another re- 
markable example 1 of strongly correlated material display- 
ing both unconventional superconducting and charge-density 
wave (CDW) phases. It is a layered compound with interca- 
lated chain and two-leg ladder subunits. Under doping and 
high hydrostatic pressure, the ladder layers exhibit supercon- 
ductivity. 2 For x — 0, the SCCO is intrinsically hole-doped 
and showed a charge modulation within the chain layers, 3 
which later was attributed to a buckling of the chains. 4 On the 
other hand, another set of experiments revealed evidence of 
commensurate charge modulations, with periods A of 3- and 
5-lattice spacings along the ladders, which bear their origin in 
the interactions between electrons. 5-7 Surprisingly, these ex- 
periments did not find modulations with period A=4 in con- 
trast to theoretical results on isolated ladders, 8 - 9 hence sug- 
gesting that the frustrated nature of the zig-zag inter-ladder 
coupling, forming a trellis lattice (see Fig. QJ), has a crucial 
role. 



In this Rapid Communication, we carry out a systematic 
investigation of the 2D charge ordering in models of cou- 
pled chains of hard-core bosons (HCB) and electronic ladders. 
We use a combination of the density matrix renormalization 
group 10 (DMRG) technique to solve the quasi-lD subsystems 
and a mean-field (MF) treatment of the screened Coulomb re- 
pulsion between them. Even though our main focus is the 
physics of doped spin ladders, the simpler hard-core boson 
model is also of great interest since (i) it displays similar 
density fluctuations and (ii) it is closely related to the ladder 
model in the limit of very large exchange rung coupling. 11 
Results from these general models are lastly compared with 
other possible interpretations of the experiments. 

Models and their ID charge fluctuations - We start by re- 
calling the behavior of charge fluctuations in an isolated doped 




FIG. 1: (a) Sketch of coupled ladders in a layer of the SCCO com- 
pound, (b) Proposed CDW modulations for A = 3, 4, 5 (as the elec- 
tronic density is identical on both leg, only one is represented). 



two-leg ladder that we describe by the isotropic t-J model: 



where V are Gutzwiller projectors, c, )CT the electron annihila- 
tion operators with spin a at site i, the density operator and 
Si the spin operator. The antiferromagnetic spin exchange J 
is set to 0.35t, a typical value for cuprates. The model has a 
spin gap that survives to finite doping S, 12 beyond S = 1/3. Its 
low-energy physics is described by a single charge mode 4>(x) 
governed by a Luttinger Liquid (LL) effective Hamiltonian 13 : 



uK(ttU(x)) 2 + — (V<Xx)) 2 
K 



(1) 



with II(x) the canonically conjugated field. Due to the pres- 
ence of a spin gap, the leading charge fluctuations are not the 
usual 2k p ones but the Akp ones. In the following, we will 
denote this wave-vector q = 27m = ikp, with n the elec- 
tronic density (S = 1 — n = 1/A). The density fluctuations 
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have a power-law decay (n(x)n(O)) c ~ cos^x)^ 2 -^ gov- 
erned by the Luttinger parameter K. Depending on K, there 
is a competition between a 4/c^-CDW (K < 1/2) and a d- 
wave superconducting (A" > 1/2) phase: the phase diagram 
has been obtained numerically. 814 For the commensurate fill- 
ings n = 3/4 and n = 1/2, two different kinds of ordered 
phases appear. 8,9 In particular, strong interactions yield a com- 
mensurate CDW for A = 4 but no signatures of translational 
symmetry breaking was observed 9 for A = 3 or 5. 

Furthermore, some insight into the 2D CDW locking will 
be obtained using the ID following HCB t-V model: 



H = -t^y>lb i+1 +hjc.] 



with bi the hard-core boson annihilation operator and V the 
nearest-neighbor repulsion. The choice is motivated by the 
analogy between hole pairs and bosons. 11 For isotropic lad- 
der, a more rigorous mapping would also involve additional 
bosons for triplets^ However, our bosonic model has the 
same effective Hamiltonian ([l} as the t-J ladder and its lead- 
ing density fluctuations are 2kp ones with exponent 2K. In 
terms of the boson density n = 1/A, we have exactly the same 
wave-vector q = 2irn = 2kp as for the ladders. In spite of this 
direct analogy, the behavior of K is different in this model. 13 
Mean-field treatment and CDW patterns - We now turn 
to the effect of coupling these quasi- ID charge fluctuations 
via screened Coulomb repulsion (we neglect other couplings, 
such as particle or pair hoppings). According to the lattice 
structure of SCCO and recent ab-initio calculations, 1 - we con- 
sider that the main interactions between adjacent ladders are 
along the grey bonds of Fig. QJ, the magnitude of which is 
denoted by V± . Using a MF approximation, these couplings 
boil down to a chemical potential term —fiiiii with 
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The local density (ny) is then determined self-consistently. 
For A = 3 and 5, the patterns which minimize the energy 
are naturally site-centered CDW shifted by tt (see Fig. QJ) 
corresponding to 

(rij) = n[l + p q cos(q(j ~ tp)) + p 2q cos(2qj)] (3) 

with <p = 0. Notice that there actually is only one har- 
monic for A = 3 and two for A = 5 with the particular- 
ity that 2q = q/2 (hence 2q = 8kp = 2kp in ladders). 
In the limit of small V±, the MF scheme is actually equiv- 
alent to static RPA as the condition for the order to develop 
is 1 = zV± r y(q)xiD{<l), with 7(g) a geometrical prefactor, z 
the coordination number and Xid(q) the static charge suscep- 
tibility of the isolated quasi- ID system. Classical configura- 
tions with one particle each A site, expected in the large V± 
limit, are also described by Eq. Q so that the MF scheme in- 
terpolates between the perturbative (RPA) regime to the non- 
perturbative (classical) one. In both repulsive models under 
study, for which K < 1, Xid(q) diverges 13 and the order 
builds up as soon as V± is branched. Focusing now on the 
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FIG. 2: (Color online) HCB model: (a) Finite size scaling of p q for 
A = 3 and V = 1. (b) results with PBC showing the large Vj_ classi- 
cal limit, (c-d) Extrapolated results with OBC showing the opening 



classical limit when A = 4, the frustration of the zig-zag cou- 
plings makes it impossible to accommodate two chains with 
one particle every four sites, as it happens for the square lat- 
tice. Striped configurations could be realized but we actually 
found that the most stable pattern is again two shifted CDW 21 
of the type (fJJ (see Fig.QJi) which can yet have different am- 
plitudes, contrary to the A = 3, 5 cases. The CDW (i) is bond- 
centered (ip — 1/2) while the other (ii) is site-centered. Here 
again, the RPA argument gives that the order develops at any 
finite V±. Numerical calculations are performed taking one 
finite quasi- ID system with open boundary conditions and 
solved by DMRG, keeping 800 states^, embedded between 
two infinite ones displaying the pattern ([3). When A = 4, the 
two patterns (i) and (ii) are solved simultaneously. In the self- 
consistent procedure, the p q 's are sampled in the bulk and the 
convergence criteria is to have a relative energy error smaller 
than 10~ 5 , for the which the density patterns are converged. 

HCB model - We start the discussion with the HCB model. 
As it can be solved efficiently, we use both periodic (PBC) 
and open boundary conditions (OBC) to check the finite size 
effects. For PBC, a finite-size chain does not display an order 
below a certain value V£(L). This can be understood in terms 
of a correlation length £ larger than L, making the scaling of 
the order parameter with L difficult with PBC. With OBC, 
there always is some charge fluctuations at the edges which 
decay toward the bulk value, the so-called Friedel oscillations. 
If the field <f> is not pined, these oscillations are well fitted by 
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FIG. 3: (Color online) Scaling of the opening of the gap for A = 3 
(a) and A = 5 (b) in the HCB model with V = 1 (extrapolated data). 



the following ansatz based on LL theory and symmetries: 



n(x) 



A 



cosq(x - (L+ l)/2) 
[sin(7ra;/(L + 1))] 



K 



(4) 



which gives an access to the Luttinger parameter K (the decay 
exponent of the 2q harmonic is 4K so that its contribution is 
negligible). When <f> gets pined, the decay of the oscillations 
are exponentially suppressed, as e~ x ^. Using OBC, we start 
the MF iterations from the V± = density profile that dis- 
plays Friedel oscillations and, once convergence is reached, 
we scale the bulk value of p q using the following ansatz which 
interpolates between the two regimes: 



p q (L) = p™ + Be- L /t/L<* 



(5) 



According to the results gathered in Fig. [2^, such an ansatz 
looks reasonable. In Fig.|2|b-d), one finds the opening of the 
order through the behavior of its Fourier components p q for 
A = 3,4,5 both in the non-interacting (V = 0) and interacting 
(V = 1) regimes. As predicted by RPA, the g-component 
opens first while the harmonic at 2q develops above V± > 0.5. 
Remarkably, p2 q hardly has finite-size effects. Notice that, for 
A = 4, the main Fourier components of the two inequivalent 
chains (i) and (ii) have the same magnitude as one may infer in 
the perturbative regime of the continuum limit. In the large V± 
limit, the Fourier components saturate close to their classical 
expectations for A = 3 and 4(i). For A = 5, the nearest- 
neighbor repulsion is not sufficient to pin the bosons every five 
sites. We found that such a classical pattern can be realized by 
taking into account longer ranged interactions within the MF 
approximation (data not shown). In all cases, we conclude that 
the patterns of Fig.QJ, similar to a Wigner crystal of bosons, 
can be realized with a sizable order for sufficiently large yet 
short range interactions. 

Although LL theory cannot give a quantitative prediction 
for the magnitude of p q , it provides the scaling of the gap 
associated with the ordering. The one boson gap is defined as 

A(L) = E {N + 1) + E Q {N - 1) - 2E {N) 

with Eq (N) the ground-state energy with N bosons. All ener- 
gies are computed with the same external /ij obtained from the 




FIG. 4: (Color online) tj model: (a) finite-size extrapolations of p q 
for A = 3. (b) Opening of the Fourier component with the inter- 
ladder coupling V± (extrapolated data for A = 3 and 5). The A = 4 
data are shown for L = 60 and cannot be extrapolated (see text). 



converged density pattern. In other words, we assume that the 
surrounding charge distribution of the neighboring chains is 
not perturbed by removing/adding one particle. Gaps are then 
extrapolated using A(L) = A°° + Be~ L ^/L. At the level 
of LL theory, the MF coupling yields a sine-Gordon term 17 
which generically reads 



g / dx cos f3<t>(x) . 



(6) 



Such a perturbation is relevant when K < K c = 8/(3 2 . For 
the main harmonic q, one has g oc V±p q and (3 = 2. Thus, 
we have K c — 2 so that when p q ^ 0, one is deep into the 
massive phase of (0. In this regime, the gap scales 23 as A oc 
g v with v = K C /2(K C - K) = 1/(2 - K). In Fig. [3^, we 
report the scaling of the gap for V = 1 and A = 3 (for which 
there is no extra harmonic), suggesting an exponent v ~ 0.84. 
A perturbative estimate of K is [1 + V/ (irt sin(7rn))] _1/2 ~ 
0.855 which gives v = 0.87. A non-perturbative estimate of 
K is deduced from fitting the density pattern using Eq. we 
find K ~ 0.78 and v ~ 0.82. The agreement with the sine- 
Gordon predictions is consequently very good up to relatively 
large V± (for these particular set of parameters). Beyond this 
regime, one recovers v = 1 which is expected when the A = 3 
pattern saturates: DMRG results nicely interpolate between 
the two regimes. In Fig.[3]3, similar results are found for A = 
5, excepting that the large V± exponent is much smaller ( v ~ 
0.28) because there is more room in the MF potential for the 
extra particle. The large V± is thus strongly dependent on the 
density and the range of the interactions. 
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t-J ladders - We now turn to the more realistic case of 
doped ladders. First, we recall that, in an isolated ladder, 
Umklapp processes can bring a sine-Gordon term (O with 
f3 = 8 when A = 4. Even though the corresponding crit- 
ical value K c = 1/8 is very small, signatures^ of this ID 
crystal has been found for small enough J / 1. From Eq. |@), 
we get the following values of K\ for J/t = 0.35 on a sys- 
tem with L = 120: K 3 ~ 0.55, K 4 ~ 0.42 and K 5 ~ 0.58. 
We notice a few qualitative differences with the HCB model. 
Firstly, K is smaller for A = 4 suggesting that this commen- 
surability will lead to the larger order in the continuum limit. 
Secondly, the coordination number is now z = 2 instead of 
z = 4 for HCB. This roughly induces a factor two in the V± 
required to stabilize a sizable order. In Fig. |4^, we give the 
finite size scaling of the Fourier amplitudes in the A = 3 case 
using ([5]). As K ~ 0.5, finite-size extrapolations are much 
harder. For A = 3 and 5, the fits give ~ 0.01, 0.02 which 
is within numerical accuracy. For A = 4, there is a finite size 
effect which makes the extrapolation impossible: for the (ii) 
chain, which pattern is ABCBABC. . . , we must take a ladder 
with one more rung to preserve the reflection symmetry about 
the center. Contrary to the HCB model (for which there is 
no issue with the finite size scaling), the minima of the den- 
sity are gradually shifted by one site between the edges and 
the middle. The consequence is that the pattern (O becomes 
frustrating at the edges, lowering the overall order on small 
systems. We observe in Fig.|4j3 that the order increases slower 
for A = 3 than for A = 4, 5. This is qualitatively different 
from the HCB model with V = 1. We may attribute this dif- 
ference to the effective behavior of pairs of holes in the ladder. 
In fact, hole pairs in isotropic doped ladders can be viewed as 
a hole pair resonating with a singlet on a plaquette. The cen- 
ter of mass of the pair can live on a rung or at the center of 
the plaquette so that their effective hard-core bosonic model 
has twice the number of sites. 11 We suggest that the slower 
opening of the A = 3 order is related to the extension of the 
hole pairs which favors overlapping at large doping. This (rel- 



ative) reduction of the density fluctuations with doping even 
leads to the absence 9 of a CDW order when n — 1/2 (A = 2). 
Another consideration is that, at the RPA level, the geomet- 
rical factor 7(g) oc cos(7r/A) favors large A. For V± > 4, 
the A = 3 curve eventually reaches the classical expectation. 
In addition to this feature, hole pairs have effective interac- 
tions at longer range 11 which show up in the low K that can 
be achieved. This should favor crystallization at small doping 
and we indeed see that the A = 5 and 4 curves are comparable. 

Conclusion - The emergence of charge ordering in a 2D ar- 
ray of two-leg ladders with a zig-zag coupling mimicking the 
crystallographic structure of the ladder planes of SCCO has 
been studied using a MF approximation to treat the screened 
Coulomb interaction. Our numerical approach provides quan- 
titative predictions for the order parameter of charge modu- 
lations with period A = 3, 4, 5. A good agreement with LL 
theory is found in the perturbative regime before reaching the 
classical limit. The analogy between the results for the two 
models suggests that such structures can be interpreted phys- 
ically as pair density waves, a localization of Cooper pairs 
with no superconducting coherence, 18 spatially organized as 
in a Wigner crystal of pairs. 6 We note that only the A = 3 and 
5 charge order have been found experimentally while our min- 
imal model still supports a sizable order for the A = 4 mod- 
ulation. An Hartree-Fock approach to a more sophisticated 
model 19 suggested a relative reduction of the A = 4 order but 
there is no clear suppression. Another proposal was to take 
into account a magnetic ring exchange 7 that lowers the A = 4 
charge order. However, according to the behavior of the pair- 
ing energy with doping and ring exchange in the t-J model, 20 
the A = 5 modulation should also be suppressed. A possi- 
ble explanation could be that the experimental set-up was not 
able to probe the peculiar type of ordering of Fig. QJ) which 
involves two inequivalent ladders. Lastly, the A = 4 pattern 
could be more sensitive to impurities, magnetic exchange, sin- 
gle particle hopping or pair tunneling between ladders leading 
to a competition with superconducting phases 
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bative regime the A = 4 staggered pattern is stabilized w.r.t. the 
stripes configuration as the interaction energy difference (per site) 
is e sta gg — e; n t oc — Pg(V2 — 1). In the classical limit where there 
is one particle each 4 site in (i), we see that the interaction energy 
is zero in the staggered pattern while it is positive for the stripes. 
The density patterns do not require a lot of states to converge, 



particularly when V± 0. For V± = 0, we checked using up to 
1600 states. 

In the non-interacting limit where K = 1, a BCS-like opening 
A oc e" 00 " 81 -/^ is expected at small V±. 



